Reentrant vortex lattice transformation in four-fold symmetric superconductors 
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The physics behind the rhombic— >square^rhombic flux line lattice transformation in increasing 
fields is clarified on the basis of Eilenberger theory. We demonstrate that this reentrance ob- 
served in LuNi2B2C is due to intrinsic competition between superconducting gap and Fermi surface 
anisotropies. The calculations reproduce not only it but also predict yet not found lock-in transition 
to a square lattice with different orientation in higher field. In view of physical origin given, this 
sequence of transitions is rather generic to occur in four-fold symmetric superconductors. 

PACS numbers: 74.60.Ec, 74.60.-w, 74.70.Dd 
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The morphology of equilibrium flux line lattice (FLL) 
in type-II superconductors, its symmetry and orientation 
relative to the crystallographic axes, is determined by 
microscopic electronic properties, particularly by Fermi 
surface topology and superconducting pair symmetry. 
In general, the shape and orientation of FLL cell will 
change with field direction trying to adjust the electronic 
anisotropy of underlying crystal M . In spite of a long re- 
search history which dates back to the pioneering small 
angle neutron scattering experiment by Cribier et al. 
on Nb, problems associated with FLL are still lively dis- 
cussed even on seemingly simplest elemental metal Nb 
||. FLL symmetry transformation under varying ap- 
plied field is one of the topics which attracted much of 
the attention of vortex physics community recently. The 
effect became known in the early 70's when the low field 
rhombic-to-square FLL transition has been observed first 
in PbTl |Q. Renewed interest in this phenomena came, 
after similar transformations has been detected in a num- 
ber of superconductors: (Re)Ni2B2C (Re=Lu,Y,Er,Tm) 
VsSi H and high-T c cuprate Lai.83Sro.i7Cu04 
All of them appear to be four-fold symmetric crys- 
tals: cubic or tetragonal. 

Recently, Eskildsen et al. jl0| discovered a remark- 
able reentrant transition of the rhombic FLL symmetry 
for H || c in LuNi 2 B 2 C (if c2 , o ~ 9 T). Upon increas- 
ing if, the rhombic lattice changes into square one and 
then backs again at a higher field. With increasing tem- 
perature the rhombus-square boundary bends away and 
never crosses if c2 line. Combining its result with other 
experiments such as STM §, 0, ^SR jjJJ and Bitter 
decoration |l| , square FLL region is confined in a small 
area just above if c i and up to ~ 2-3 T on (if, T) plane. 
High field square-to-rhombus transition is detected in 
TmNi2B2C as well, though this is a magnetic member 
of borocarbide family. It is quite interesting to remember 
Nb case. The FLL in this cubic system for H || (001) 
exhibits square-to- rhombic transition as T increases [O . 
Although the definite phase diagram is not established 
yet, temperature dependence of FLL symmetry alone 
suggests an isolated region of stable square lattice in if 
vs. T. 




It was realized early [19] that a certain four-fold 
anisotropy in plane perpendicular to the applied field, 
such as the Fermi velocity vp, can drive low field 
rhombus-to-square transformation. In the in-plane 
anisotropy modeled by Vp{9) = vp(j){l + (3 cos 40) (9 
polar angle), square FLL D v with nearest neighbors ori- 
ented along the velocity minimum will be stabilized in low 
fields if anisotropy degree /3 is large enough. The same is 
true for the four-fold gap anisotropy alone |A(r,#)| 2 = 
A(r,-|)| 2 (l — acos46>), when square FLL D g (nearest 
neighbors along gap minimum) tends to be stabilized. 
These are schematically shown in Fig.l. The nontriv- 
ial question is why square FLL is unstable at high fields 
and what the actual phase diagram looks like? Contrary 
to the answer given by Gurevich and Kogan JT(| who 
consider it due to thermal fluctuation near if C 2, here we 
investigate two indispensable anisotropic effects on the 
same footing. The interplay of gap and Fermi surface 
anisotropy indeed can give rise to the reentrant FLL tran- 
sition and further square lock-in transition in a higher 
field. 

There is firm evidence in LUM2B2C and YM2B2C, 
which we regard as essentially the same systems, to 
show the existence of both anisotropies. As for the gap 
anisotropy, various bulk measurements, including the in- 
dependent linear specific heat coefficient 7 (if) (Volovik 
effect) fL7| , the thermal conductivity JTsJ, and angle- 
integrate photo-emission |fl9|| , all show a substantial gap 
anisotropy. More recently the angle-resolved thermal 
conductivity [^0) under if exhibits that the gap is van- 
ishingly smaller in (100) than in (110). This is consis- 
tent with the spatial extension of the zero-bias peak ob- 
served by STM |^l| . As for the Fermi velocity anisotropy, 
when interpreted through simple four- fold harmonic vari- 
ation, band calculation [p2| told us that the Fermi ve- 
locity is larger in (100) than in (110), or the angle re- 
solved density of states (DOS) N(9) ~ l/vp(9) is smaller 
along (100) compared to (110). Note that the sense of 
two anisotropies mentioned is what we would naively ex- 
pected because the larger gap A(9 = 0) should develop 
in larger N(9 = 0). The built-in tendency to stabilize 
two different orientation □ „ and D g of square lattice ul- 
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FIG. 1: Rhombus and square FLL cell orientations relative to 
the crystal. The vortex center is shown by •. a > 0, /3 > 0. 



FIG. 2: FLL phase diagram at T/T c = 0.5. a) /3-dependence 
(a — 0), o shows where □„ is stable, b) a-dependence (f3 = 
0.3). Broken lines are guide for the eye. 



timately leads to the rich vortex phase diagram, including 
the reentrant as we will see shortly. 

Apart from the limiting cases: London model at H <C 
H c2 (T) and Ginzburg-Landau (GL) model at T w T c , 
there is no handy and convenient approximate scheme 
to describe microscopic vortex properties deep in (H, T) 
plane. Therefore we resort to quasi-classical Eilenberger 
equations J23| valid for kp^a 3> 1 (kp: the Fermi wave 
number and £o : the coherence length) a condition met in 
most of superconductors. Eilenberger equations read as 
(71=1) 

(2uj + v F (0) • II) f(w, r, 9) = 2A(r, 0)g(u, r, 9), (1) 

(2w - v F (9) ■ IT) /t(w, r, 9) = 2A*(r, 0)g(u, r, 9). (2) 

Here II = V + (2m /^^) A is gauge invariant gra- 
dient, A is vector-potential and $o is flux quantum. 
uj = nT(2n + 1) with integer n is Matsubara fre- 
quency. Normalization condition for Green's function 
g 2 + ff> = 1. The pairing interaction is assumed sep- 
arable V(9,9') = V 0(9) 0(9') so that gap function is 
A(r,9) — *f?(r)(j)(8). We consider two-dimensional case 
with cylindrical Fermi surface. Four-fold models for 
Fermi velocity vf{9) = + /3cos4#)/yl — (3 2 and 

gap anisotropy (f)(9) 2 — (j>(j) 2 (l — acos4#) have been 
adopted. Here 9 is polar angle relative to (100) axis. 
Constant 0(f) 2 = 1/(1 + (1 - a/I - (3 2 )a/(3) is cho- 
sen to assure the same T c and DOS Nq for any value 
of anisotropy parameters a and (3. We are most inter- 
ested in case when a and (5 are of the same sign. Then 
positions of gap and velocity minima are 45° rotated to 
each other (see Fig.l) and we have a competing effect. 
The self-consistent equations for the gap function \&(r) 
and vector-potential A are 
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V x V x A = -^-N Tv F Im (gu) ■ (4) 



Here, v(9) = vf(9)/vf and u — (cosfl, s\w9) is unit vec- 
tor along the Fermi velocity vp = vp(9)u. For average 
over Fermi surface, (...) = (l/27r) j . . . d9. Extra fac- 
tor l/v(9) in averages came from angle resolved DOS 
N(9) = N /v(9) on Fermi surface. 

We have performed extensive numerical computations 
by the so-called explosion method (see Ref. ^4] for details) 
for various values of anisotropy parameters a and /3 in a 
high GL parameter case k=100. The self-consistent solu- 
tion yields a complete set of the physical quantities: the 
spatial profiles of the order parameter ^(r) and the mag- 
netic field H(r), and the local density of states around a 
vortex core. The free energy density is given by 



F = ^ - 2nTN ^(^m^±in) . (5) 
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Here, a — (B/$o) f ce u a dr. Free energy should be min- 
imized with respect to the FLL symmetry and its ori- 
entation relative to the crystallographic axes. Numerics 
is backed up by analytical calculations. Namely we also 
solve these analytically at the two limiting cases; (1) near 
H c2 (T) and (2) H < H c2 (T) to gain physical insights. 
For analytical results we considered FLL cell shaped as 
rhombus with apex angle in interval [60°, 90°]. Two dif- 
ferent orientations of rhombus cell are compared: rhom- 
bus diagonals along gap minimum (velocity maximum) 
and along velocity minimum (gap maximum) (see Fig. 
|l|). Since the numerical computation is very demand- 
ing and time consuming, we limit ourselves to the four 
configurations (apex angles 60° and 90° for each of two 
orientations) as shown in Fig.l. 

In Fig. ||(a) we show the phase diagram for t = T/T c = 
0.5 and a — 0, where the square lattice □„ becomes sta- 
ble if the Fermi velocity anisotropy /3 exceeds a certain 
value. It shows also that the triangular lattice is always 
stabilized at lower and higher H . Thus even without the 
gap anisotropy (a = 0) the reentrant transformation can 
be induced. This can be understood by looking at the 
self-consistency Eq. (|3|). It is seen that factor \/v(8) 
plays the role of gap anisotropy even in isotropic pairing 
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3: Field dependence of free energy for A v , A g and O g 
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case (0(0) = 1). 

The gap anisotropy further induces a rich variety of 
the phase diagram. In Fig. ||(b) we show it when the 
moderate velocity anisotropy (3 = 0.3 is taken. We see 
that at the large a cases the triangular lattice directly 
changes into the square lattice D g . This is a similar 
situation to the <i-wave case with the isotropic Fermi 
velocity |p3. The preferred orientation perfectly coin- 
cides with the nodal direction, namely, nearest neigh- 
bors of FLL are along (100) (along (110)) in d xy (d x 2_ y 2) 
pairing. As a decreases, this d g lattice region shrinks. 
At lower fields yet another transformation from the tri- 
angular to square rotated by 45° relative to n g , 
emerges. The most complicated case is at moderate val- 
ues of a and (3, leading to successive lattice transforma- 
tion A — >• □„ A ^ D 9 as iJ increases. In a weaker a 
case the last lock- in transition is absent. 

As a typical example we investigate the case a — (3 — 
0.3 in details at lower temperature. In Fig. ||the four free 
energy curves are shown as a function of h = H/H C 2, o 
at t = 0.25. As h increases the above mentioned suc- 
cessive transformations are clearly illustrated as several 
crossings. Figure ^ displays the resulting phase diagram 
in (H, T) plane. It is seen from this that (1) the □„ 
region is confined to lower H and T; (2) its boundary 
bends away from H C 2\ (3) the high field region is occu- 
pied by the d g lattice; (4) along H C 2 line this terminates 
at t = 0.56, below which the rhombic lattice becomes 
stable. 

Let us discuss the physical origin of this intricate phase 
diagram in connection with the observation in borocar- 
bides. Basically, in order to induce narrowly limited □„ 
region at low H and T we need the competing effects, 
each coming from the gap and velocity anisotropies. The 
gap anisotropy with a > prefers the lattice D g with 
nearest neighbors along (100) direction, while the veloc- 
ity anisotropy with j3 > tends to favor 45° rotated 
square lattice Thus these anisotropies compete each 
other. The a(3 < case does not cause such a frustration 
in FLL. 

The interplay of two anisotropies can be understood 




FIG. 4: Phase diagram of FLL for a = (3 — 0.3. In inset, 
apex angle variation along H C 2 line. The solid (broken) line 
is for orientation "v" ("g")> see Fig. EL 



by looking at the free energy density at lower H: 
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A is penetration depth, and ^Sq(T) is zero field gap. The 
extra four-fold anisotropy due to non-local correction ap- 
pears through the parameter d — t xxxx + t yyyy — 6t xxyy in 
the above. At lower T, tiji m contains factor v 3 (0) / cb 3 (0) 
which strengthens the tendency towards the velocity 
anisotropy ((3 is effectively increased by a when a/3 > 0). 
Near T c this factor becomes <j) 2 (8)v 3 (8) which weakens 
the combined anisotropy effect by canceling each other. 
Therefore the two anisotropies plays different role, de- 
pending on T, giving rise to the bent transition curve 
in the (H, T) plane. For a fixed temperature as H in- 
creases from H&, FLL starts with the regular triangle 
lattice because the electromagnetic interaction between 
far apart vortices is isotropic, yielding the closed packed 
hexagonal symmetry just above H c \. As vortices ap- 
proach each other, proliferating anisotropy in current dis- 
tribution will squeeze rhombus FLL cell toward square 
shape. Being dominant at low T region, anisotropy of 
Fermi velocity will stabilize square lattice □„ as soon 
as vortices come close to each other with increasing H. 
The observed orientation of D v with nearest neighbors 
along (110) is indeed expected in LUNI2B2C since the 
band structure calculation indicates that Fermi velocity 
is larger in (100) than (110) when mapped into our four- 
fold model p2| . The same result is obtained within non- 
local London model with Fermi velocity anisotropy alone 
p5[ . Upon further increasing H, condensation energy 
gradually takes over the major role in determining the in- 
teraction between vortices. Physically it is due to the ki- 
netic energy cost of quasi-particles localized around core 
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(see Ref. g4| for details). In high fields, at least in high-K 
superconductors, the anisotropy in vortex- vortex interac- 
tion is exclusively due to the vortex core anisotropy giv- 
ing rise to the transition from the low field □„ to high field 
□ g via intermediate rhombic lattice (□„ — > A — > D g ). 

Limit of H « i? c2 (T) allows the analytical solution for 
free energy 
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Here, /, /' and \I/ are solutions of the linearized Eilen- 
berger equation, h s is magnetic field due to supercur- 
rents, B = H + h s is magnetic induction. In the isotropic 
case a = (3 = and near T c it reduces to the standard 
Abrikosov expression T — (2 k 2 — 1)0 a- This expression 
for free energy is valid all along the H C 2(T) line. The 
result is illustrated in Fig. 4 as the inset showing that 
the apex angle of the rhombic lattice continuous changes 
about the regular triangular lattice with 60° . The diago- 
nal of the rhombic lattice rotated by 45° at t = 0.64 from 
(100) to (110), i.e. from to A 9 . At a higher field A s 
lattice is locked in D g . This lock-in point t = 0.56 is 
rather nicely connected to the points determined numer- 
ically as shown in the main panel in Fig. 4. This implies 
that our numerical results, which examine the limited 
number of lattice configurations (4 types), yield a rea- 
sonably phase diagram even taking into account general 
rhombic lattices. 

The present calculation does not aim to quantitatively 
reproduce the actual FLL phase diagram in LuNi 2 B 2 C, 
but to physically understand its possible lattice transfor- 
mation. Note that so far the observed square lattices in 
LuNi 2 B 2 C by several methods are all not n g . The D g 
phase can occur at further high H if the gap anisotropy 
is large enough. Since this is the case for LuNi 2 B 2 C as 
mentioned before, there is a good chance to observe it. 
According to various thermodynamic and transport ex- 
periments the gap anisotropy is extremely large, the re- 
cent directional dependent measurements shows the node 
situated at (100) direction. This unambiguously tells us 
that LuNi 2 B 2 C belongs to the most interesting compet- 
ing case (a > and (3 > 0) and also gratifying enough 
that it belongs to the naively expected case where maxi- 
mum energy gap is oriented to the maximum DOS direc- 
tion. So far the investigated field for the vortex structure 
is limited below 4T (compare i? c2 , ~ 9 T). 

As for the other materials where the square lattice is 
found such as Nb @, V 3 Si j§ and Sr 2 Ru0 4 @ we can 
expect the reentrant transition from the square to rhom- 
bic lattice and further lock-in transition to the square 
lattice with different orientation if the gap anisotropy is 



strong enough. The recent finding of rhombus-to-square 
transition at lower H in La1.83Sro.17Cu.O4 |J deserves 
a special attention because in spite of d x 2_ y 2 symmetry 
they discovered not D g . This nontrivial observation 
is indeed expected by our calculation p7| . 

In conclusion we have shown that the vortex lattice 
morphology is deeply connected to the underlying mi- 
croscopic electronic structure. Specifically it is seen that 
the reentrant transition from the square to rhombic lat- 
tice in LuNi 2 B 2 C can be well understood as arising from 
the two competing superconducting gap and Fermi sur- 
face anisotropies both of which are documented to exist 
experimentally. It is demonstrated by solving the quasi- 
classical Eilcnbcrgcr equations numerically and analyti- 
cally. We have shown another yet not found FLL trans- 
formation to the differently oriented square lattice in a 
higher field. In view of physical origin given, this kind of 
re-entrance and high field square lattice are rather generic 
to occur. Thus we expect a similar successive transition 
in type-II superconductors with four-fold symmetry at 
least with large k. 

We thank Y. Matsuda, N. Nishida, K. Takanaka and 
E. M. Forgan for valuable information. 
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